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Astrophysical interpretation of the medium scale clustering in the ultrahigh energy sky 
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We compare the clustering properties of the combined dataset of ultra-high energy cosmic rays 
events, reported by the AGASA, HiRes, Yakutsk and Sugar collaborations, with a catalogue of 
galaxies of the local universe (redshift z < 0.06). We find that the data reproduce particularly well 
the clustering properties of the nearby universe within z < 0.02. There is no statistically significant 
cross-correlation between data and structures, although intriguingly the nominal cross-correlation 
chance probability drops from O(50%) to O(10%) using the catalogue with a smaller horizon. Also, 
we discuss the impact on the robustness of the results of deflections in some galactic magnetic 
field models used in the literature. These results suggest a relevant role of magnetic fields (possibly 
extragalactic ones, too) and/or possibly some heavy nuclei fraction in the UHECRs. The importance 
of a confirmation of these hints (and of some of their implications) by Auger data is emphasized. 
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I. INTRODUCTION 



One of the keys towards the solution of the mysteri- 
ous origin of ultra-high energy cosmic rays (UHECRs) 
is the study of their anisotropy pattern. The chances 
to perform (some kind of) UHECR astronomy increase 
significantly at extremely high energy, in particular due 
to the decreasing of deflections in the galactic and pos- 
sibly extragalactic magnetic fields. Moreover, at E > 
(4 — 5) x 10 19 eV the opacity of the interstellar space 
to protons drastically grows due to the kinematically al- 
lowed photo-pion production on Cosmic Microwave Back- 
ground (CMB) photonSiknown as the Greisen-Zatsepin- 
Kuzmin or GZK effect [H, HJ . A similar phenomenon at 
slightly different energies occurs for heavier primaries via 
photo-disintegration energy losses. Recently, an observa- 
tional evidence for a flux suppression consistent with the 
GZK feature has been reported by the HiRes collabo- 
ration 0]. Within reasonable astrophysical assumptions, 
these energy-losses phenomena impose a conservative up- 
per limit to the distance from which the bulk of UHECRs 
is emitted, of the order of a few hundreds Mpc at most, 
which may enhance the chances of identifying structures. 
In Ref. Jj|a forecast analysis for the Pierre Auger Obser- 
vatory [5], lf| was performed to derive the minimum statis- 
tics needed to test the "zeroth order" hypothesis that 
UHECRs trace the baryonic distribution in the universe. 
Assuming proton primaries, it was found that a few hun- 
dred events at E > 5 x 10 19 eV are necessary at Auger to 
have reasonably high chances to identify the signature. 
On the other hand, available catalogues from the exper- 
iments of the previous generation contain 0(100) events 
above E > (4 — 5) x 10 19 eV, thus motivating a search 
for possible angular patterns already in the present data 
In particular, after renormalizing the energy scales 
of the different experiments to the HiRes one at 4 x 10 19 



eV, the authors of @ found some evidence of a broad 
maximum of the two-point autocorrelation function of 
UHECR arrival directions around 25 degrees. Since the 
search was made a posteriori, the assessment of its signif- 
icance is a delicate issue and involves the determination 
of a penalty factor critically dependent on the performed 
number of trials. The previous claim of small scale clus- 
tering in the AGASA data and the associated debate on 
its significance (see e.g. d, [Io|, HH ) would suggest to take 
a cautionary attitude towards a posteriori claims. How- 
ever, a similar feature has been found in the Auger data 
alone as well, as recently reported in [T3|. We shall thus 
proceed in the following under the assumption that the 
signal is real, exploring some astrophysical implications. 

In [l3[ , the present authors already tested the qualita- 
tive interpretation of the result (as reflecting the large- 
scale structure (LSS) of UHECR sources) given in Q on 
the light of our previous map templates obtained from 
the IRAS PSCz galaxy catalogue [T3]. The observed 
data and the Monte Carlo events from the catalogue 
share several features, which are even more prominent 
if a quadratic correlation with LSS is assumed. On the 
other hand, no relevant cross-correlation has been found, 
which would be the smoking gun to test such scenarios. 
However, this is not particularly surprising: apart for 
the sake of simplicity, there is no a-priori reason to ex- 
pect that cosmic rays are 100% made of protons, that the 
effects of magnetic fields are negligible above 4 x 10 19 eV 
for the angular scales considered, and that the sources 
trace in an unbiased way the LSS. If some or all these 
assumptions are relaxed, the possibility of a consistent 
scenario emerges: At "low" energy, both clustering and 
cross-correlations in UHECRs are absent, since magnetic 
deflections and a very large energy-loss horizon destroy 
them. With sufficient statistics and at sufficiently "high" 
energy, cross-correlations should eventually emerge, both 
because magnetic deflections scale like 1/Encrgy and be- 
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FIG. 1: Skymap of the UHECR arrival directions of events 
in galactic coordinates with rescaled energy E > 4 x 10 19 eV. 
The solid line is the celestial equator. 

cause of the expected shrinking of the horizon. Before 
this stage is reached experimentally, it is likely that the 
first hint will appear in the clustering, but not in the 
cross-correlation. The reason being that the former is 
much more robust versus magnetic deflections than the 
second one, as we shall argue. 

In this paper, we extend our previous analysis in two 
ways: (i) we assume a smaller horizon, i.e. biasing the 
correlation with LSS towards closer sources; (ii) we study 
the impact of the galactic magnetic field (GMF) on the 
autocorrelation signature and on the cross-correlation 
signal. We anticipate that the data reproduce particu- 
larly well the clustering properties of the nearby universe 
within z < 0.02 and they are also quite robust with re- 
spect to deflections in galactic magnetic fields. We sum- 
marize our assumptions and techniques in Sec. |IIJ while 
devoting Sec. IIIII to present our results and attempting 
some interpretations of them. In Sec. IIVI we briefly dis- 
cuss our findings and conclude. 

II. ASSUMPTIONS AND METHODS 

A. The data 

In our analysis, we closely follow the approach reported 
in d, [TJ[ , using a similar dataset extracted from available 
publications or talks of the AGASA [HI Yakutsk 
SUGAR [TJ], and HiRes collaborations [U, [TJ|. Note 
that we rescale a priori the energies of the experiments to 
the HiRes one and consider events above E > 4 x 10 19 eV 
in this renormalized sample. This approach is applied 
hereafter in the analysis and we address the reader to @ 
for further details. In Fig. [TJwe show the points used in 
this analysis in galactic coordinates, while Fig. [2] reports 
the single and combined exposure for the various experi- 
ments as a function of the declination, in the limit of sat- 
urated acceptance and mediated over the right ascension 
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FIG. 2: Single and combined exposures for the various exper- 
iment considered: Sugar (red), Hires (green), Agasa (blue), 
Yakutsk (yellow), combined (black). 
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FIG. 3: The UHECR flux contrast map (or excess map) prop- 
erly smoothed with a Gaussian filter of 10° width. 

(see e.g. [Io[). In Fig. [3] we show the derived UHECR 
excess map (flux over average expected flux, minus one) 
properly smoothed by a gaussian filter of 10°. Such a 
choice for the width amplitude (which has only illustra- 
tive purposes) represents an acceptable compromise be- 
tween the few degrees of the experimental uncertainty on 
the arrival direction of UHECR, and the typical angular- 
length of the nearby astrophysical structures of several 
tens of degrees. Of course, the data have been properly 
weighted by the exposure. 

The smoothed map in Fig. [3] clearly shows that the 
most apparent visual feature in the data is the medium 
scale clustering, with the data clustered in few spots of 
20°-30° degrees each, that in their turn are distributed 
almost uniformly in the sky. This is of course the reason 
of the signal found in 0] with a proper statistical anal- 
ysis. The clustering seen in the southern hemisphere is 
due to the Sugar data only and has thus a weak statis- 
tical evidence. However, the clustering signal is indeed 
present and statistically significant also considering the 
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data from the northern hemisphere only Q . Indeed, hints 
of this clustering in the northern hemisphere were recog- 
nized already some years ago [2l|. Finally, lowering the 
energy threshold the signature disappears, excluding the 
possibility that the signal is only a systematic feature 
coming from an incorrect modeling of the exposures [|[ . 

Concerning possible interpretations of the signal, here 
we report a few general considerations, while more quan- 
titative analyses are reported in the following. The ab- 
sence of a correlation with the Galactic Plane or the ab- 
sence of excess toward the Galactic Center disfavor re- 
spectively Galactic astrophysical sources and heavy relic 
decays in our Galactic Halo as origin of these events. 
Qualitatively, extragalactic sources of astrophysical na- 
ture appear the most likely accelerators. In these scenar- 
ios, unless only a handful of sources dominate the emis- 
sion, the pattern of the arrival directions of the events 
should reflect to some extent the one of large scale struc- 
tures in the nearby universe. Actually the degree of 
clustering observed is quite pronounced. It exceeds also 
the anisotropy expected in the minimal case of proton 
primaries (with a GZK horizon z ~ 0.06), which is in 
marginal agreement with the data (see [l3| and section 
IIIip . Thus, the few prominent structures visible naturally 
suggest either a scenario where the UHE sky is dominated 
by few nearby powerful sources or one where UHECRs 
are produced by a relatively larger number of sources sig- 
nificantly biased with overdensities in the local universe 
(within z ~ 0.02). Both scenarios require an important 
role of magnetic fields, either galactic or extragalactic, 
to accommodate deflections of the order >10°. In the 
former case such deflections are necessary to explain the 
large smearing of the point sources emission and in the 
latter to justify the lack of a significant correspondence 
between the data and the nearby galaxy clusters. 

B. The models 

The previous discussion motivates an extension of the 
previous analysis reported in [l3| along two directions: 

(J) assuming a smaller horizon, i.e. biasing the corre- 
lation with LSS towards closer sources; 

(II) studying the impact of the GMF on the autocorre- 
lation signature and on the cross-correlation signal. 

Both extensions should be regarded as "first order" re- 
finements of our previous study. The point (//) does not 
need much justification: it is important to establish the 
robustness of the previous results with respect to the ef- 
fects of astrophysical magnetic fields. Even if extragalac- 
tic magnetic fields may have a major role in shaping or 
preserving the UHECR anisotropies, very little is known 
about them (see [22| and [HI]). On the other hand we 
know for sure that a regular GMF exists — although we 
have only rough ideas on its magnitude and structure — 
and it is in principle relevant for UHECR deflections, 



even in the case of pure proton composition. In the fol- 
lowing, we shall then consider how our results change 
when data are corrected for the effects of a few GMF 
models available in the literature. In particular, we shall 
use the three models HMR, TT, and PS employed in 
[2^ |. which we refer to for details. To account for the 
deflections in the GMF in our analyses we shall follow 
the back-tracking technique described in [25|, [26|. The 
technique consists in mapping the arrival CRs directions 
on the Earth backward outside the GMF to obtain a 
map of the GMF deflections. We then apply the map- 
ping to the extragalactic expected CR map and correct 
it for the GMF displacements. The extragalactic CRs 
map F(E cu t, Cl) expected at an energy greater than -E cu t 
at the direction Cl is obtained as described in The 
map is then convolved with the GMF deflections to have 
F(E cut , Cl(Cl')) where Cl(Cl') is the mapping produced by 
the back-tracking technique. This method is fully suited 
for the cases in which energy losses along the particle 
track are negligible and when the particle energy is large 
enough to exclude loops and/or trapped regions during 
the propagation. Both these conditions are satisfied for 
the UHECRs and for the GMFs we considered. Also note 
that an isotropic sky remains isotropic under the GMF 
transformation, in agreement with the ex pec tation from 
the Liouville theorem. We refer to [24l [25l.l26| for details. 

For simplicity we only consider the mapping produced 
for a fixed rigidity corresponding to the energy E cut 
(with the choice -E cu t=40 EeV in the present case). This 
should be a reasonable approximation, equivalent to re- 
place the steep (oc E~ 3 ) UHECRs spectrum above E cnt 
with a delta-function at -E C ut- Beside the steepness of the 
UHECR spectrum, a further motivation for this approx- 
imation is that we are not considering the shift of single 
objects but of an overall map, already smoothed at a scale 
of order ~ 5° (we are not interested to the very small 
scales, indeed). Finally, we shall show in the next section 
that the auto-correlation analysis is quite insensitive to 
the details of the GMFs or the assumed rigidity as long 
as the magnetic deflections effects remains moderate, so 
that the approximation is also justified a posteriori, at 
least for autocorrelations studies. The effects are poten- 
tially larger for cross-correlation analyses, which however 
are already less robust for other reasons. 

A further possible problem is given by the fact that the 
mask region present in the catalogue and excluded from 
the analysis is distorted by the effect of the GMF, so that 
in principle one should exclude, case by case, the regions 
which the mask is mapped into by the GMF. We ne- 
glect this effect assuming the the mask is approximately 
mapped in itself by the GMF transformation. This is a 
quite good approximation for the region near the galactic 
plane while it is not satisfied by the two narrow stripes. 
However, the the stripes amount to about 10% of the to- 
tal mask and only roughly 2% of the whole sky, which is 
a very small bias for our purposes in this work. 

From a qualitative point of view, the point (I) is rea- 
sonable, too. In many scenarios, only relatively nearby 
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FIG. 4: Top row: Excess maps of PH (left) and SH models (right) in galactic coordinates. The grey contour bounds the blind 
region of PSCz catalogue. Bottom row: Galactic excess maps of the PH (left) and SH models (rigth) taking into account the 
Galactic magnetic field correction for a specific model (HMR in [24[) and assuming proton primaries (Z = 1). 



sources (if any) may be identifiable in cosmic ray maps. 
There are several plausible reasons for that. Even in ab- 
sence of magnetic fields, an heavier composition implies 
a different energy-loss horizon for UHECRs [2?], HH . For 
the energy threshold considered here (E > 4 x 10 19 eV), 
this is smaller than the proton one. In presence of ex- 
tragalactic magnetic fields, the propagation of a UHECR 
may greatly differ from a straight line and in principle 
may even happen in a diffusive regime [H, [3(| HH, H3] . 
Although it is unlikely that the propagation is truly dif- 
fusive, Gpc-scale pathlcngths for protons injected within 
a few hundreds Mpc may be common even above 10 EcV 
[33l ] . A non- negligible role of magnetic fields would have 
two consequences: for a given energy-loss mechanism, it 
is clear that the true horizon may be significantly shorter 
than the expected one. Thus, UHECRs above a given E t h 
may be largely collected within a region smaller than the 
linear energy-loss horizon. More important, apart for en- 
ergy losses, the longer the propagation time, the smaller 
the chance that intrinsic anisotropies may survive (in 
some form). Finally, since UHECR source likely have to 
meet special accelerator requirements, it is reasonable to 
conceive a relatively rare population of sources, possibly 
strongly biased with respect to LSS. 

However, how to implement in practice point (J) is ad- 
mittedly not model independent. One possibility may be 
to cut arbitrarily a LSS catalogue to some redshift z cu t, 
and consider only correlations with structures within this 



distance, assuming for the rest that UHECRs are un- 
biased tracers of LSS (i.e. neglecting otherwise energy 
loss effects). Another possibility is to create anisotropy 
map templates of specific scenarios for UHECR compo- 
sition, sources, and extragalactic magnetic fields, com- 
paring them with the observed configurations of data in 
order to infer the best model. Although this will be the 
way to proceed when high statistics will be achieved, at 
the moment it could just dilute the basic consequence 
of our assumption (/) under a large number of unknown 
parameters. To keep some physical-inspired input in a 
toy model, we shall compare the distribution of data as 
in Fig. [1] with the LSS maps obtained by convolution 
of the PSCz catalogue with an energy-loss window func- 
tion corresponding to protons twice more energetic, i.e. 
E = 8 x 10 19 eV, implying an effective horizon z ~ 0.02 
[U. We shall denote this scenario as "small horizon" 
(SH), as opposed to the "proton horizon" (PH) as treated 
in jl 31 ] and corresponding to the minimal assumption of 
protons primaries with E > 4 x 10 19 cV propagating in 
a negligible EGMF (usual GZK horizon z ~ 0.06). In 
the top two panels of Fig. [4] we report the PH and SH 
maps. The smoothing is variable and it is related to the 
adaptive smoothing applied to the PSCz catalogue to 
minimize the effect of the shot noise. We emphasize that 
this should be considered a toy model, and not a realistic 
scenario for UHECR sources or composition. However, 
our toy model may be indicative of a plausible situation 
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where at least for anisotropy searches the useful UHECR 
horizon is relatively short. 

C. Statistical tools 

For the statistical analysis wc define the (cumulative) 
autocorrelation function w as a function of the separation 
angle 6 as 

N d i-1 

«>(<*) = ££e(*-^), (i) 

i=2 j=l 

where is the step function, Nd the number of 
CRs considered and 5ij = arccos(cos pi cos pj + 
sin pi sin p 3 - cos(0i — <pj)) is the angular distance between 
the two cosmic rays i and j with coordinates (p, 4>) on 
the sphere. Analogously, one can define the correlation 
function £(<$) as 

N d N s 

£(*) = EE (*-*«)' ( 2 ) 

i— 1 a— 1 

where 8i a is the angular distance between the CR i and 
the candidate source a and N s is the number of source 
objects considered. 

We perform a large number M ~ 10 5 of Monte Carlo 
simulations of A data sampled from a distribution on 
the sky corresponding to the hypothesis H (e.g., uniform, 
LSS, etc.) and for each realization j we calculate the au- 
tocorrelation function Wj(5). The sets of random data 
match the number of data for the different experiments 
passing the cuts after rescaling, and are spatially dis- 
tributed according to the exposures of the experiments. 
The formal probability P H (8) to observe an equal or 
larger value of the autocorrelation function by chance 
is 

M 

P H {S) = -Y,^f^)-^% (3) 

i=i 

where w*(5) is the observed value for the cosmic ray 
dataset and the convention 6(0) = 1 is being used. Rel- 
atively high values of P and 1 — P indicate that the data 
are consistent with the null hypothesis being used to gen- 
erate the comparison samples, while low values of P or 
1 — P indicate that the model is inappropriate to ex- 
plain the data. That is, in the following we shall plot the 
function P(6) x [1 — P(6)], which vanishes if any of P or 
1 — P vanishes and has the theoretical maximum value of 
1/4. Thus, the higher its value is the more consistent the 
data are with the underlying hypothesis. Note also that 
by construction the values at different 5 of the function 
P(5) are not independent. 

To calculate the cross-correlation probability, we per- 
form a large number M(~ 10 5 ) Monte Carlo realization 
of N events is sampled according to the LSS probabil- 
ity distribution, and for each realization i wc calculate 



the function £^ (S). We generate analogously M ran- 
dom datasets from an uniform distribution, and calculate 
£" m ((5). We have thus M 2 independent couples of func- 
tions (i, j). The fraction of the M 2 simulations where the 
condition £ uni (<5) > £ LSS (c5) is fulfilled is the probability 

M M 

p ? oo = M2 E E e KT (*) - ts 33 ®] ■ (4) 

i=l 3 = 1 

A technical detail of the analysis is related to the pres- 
ence of the catalogue mask. This includes a zone centered 
on the galactic plane and caused by the galactic extinc- 
tion and a few, narrow stripes which were not observed 
with enough sensitivity by the IRAS satellite. These 
regions are excluded from our analysis with the use of 
the binary mask available with the PSCz catalogue it- 
self. This reduces the available sample by about 10%. 



III. RESULTS 

By repeating our analysis in [l3| following the SH 
model and without considering for the moment the effects 
of the GMF, we obtain the results shown in Fig. [5] The 
SH model seems to explain extremely well the clustering 
properties of the data, with the related P x (1 — P) curve 
almost coincident with the ideal P = 0.5 expectation. 
Not surprisingly, this can be understood after a visual 
inspection of the maps in Fig. [3] (data) and Fig. [1] (mod- 
els). While the map from protons with E > 4 x 10 19 eV 
(the PH model) is still too much isotropic with respect 
to the data, in the SH map the number of clusters and 
their distribution resemble much more the data, in that 
it leaves typical "voids" between clustered hot-spots ob- 
served. 

Our next step is to investigate the effects of the galac- 
tic magnetic field. In the two bottom panels of Fig. 2] we 
show the effective modification of LSS structures happen- 
ing for the PH and the SH scenario, assuming as example 
the HMR model in [24| , In general, besides the shifting 
of the positions of the structures, as expected the GMF 
introduces in the deflected maps also other peculiar lens- 
ing phenomena like shearing and (de)magnification [2(| . 
More quantitatively, the effects of the GMF are studied 
in the following through the modifications induced in the 
auto and cross-correlation functions. In Fig. [6] we inves- 
tigate the effect of the GMF on the signature in the au- 
tocorrelation function. Wc also show the effect of chang- 
ing the rigidity of the particles. Actually, the equations 
of motion for CRs in the GMF only depend on the pa- 
rameter C — B x Z/E where B, Z, E are respectively the 
GMF normalization, the particle atomic number (electric 
charge of the nucleus) and the particle energy. A com- 
bination of parameters that leaves unchanged C is thus 
completely degenerate from the point of view of propa- 
gation in the GMF (though, of course, not for the energy 
losses in the propagation in the extragalactic sky.) As a 
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FIG. 5: Chance probability of auto correlation taking as reference model an uniform distribution, the linear correlation model 
of [TH (LSS-PH model) and the presently considered model with a smaller horizon (LSS-SH model). 



general consideration, we see that at least for the base- 
line cases considered, the correction for the GMF does 
not destroy the pattern in P(l — P), but can improve 
or worsen it at most by a factor of a few. On the other 
hand, extreme changes in C may significantly alter the 
pattern of the function. 

Finally, in Fig. [7] we report the results of the cross cor- 
relation analysis. Note that, while in the previous case 
we were repeating the same test of Ref . Q , here we per- 
form a different test, and to assess the confidence level 
of any statistically significant signal we might find one 
should carefully evaluate the penalty factor. Unfortu- 
nately, in no case we find a statistically significant signal 
(namely, not even at the nominal level). Yet, qualita- 
tively all the SH maps show some improvement with a 
nominal P^ ~ 10% (with respect to Pf. ~ 50% for the 
PH case) , even without use of the GMF correction. This 
is understood since we have about the same number of 
clusters in the data and in the map and typically one can 
find a correspondence between the two within a radius of 
roughly 50°. A more significant signal of cross-correlation 
should eventually peak within ~ 25° (the typical size of 
the clusters) signaling a superposition of the data and 
map clusters. Once again, no significant difference arises 
when the data are corrected for the GMF. Although the 
minimum of the probability can change by up to a factor 
of a few, it does not move towards S ~ 0°, as it should 
be if the GMF were correctly shifting the hot-spots. 

An important point to stress is that while an evidence 
of cross-correlation would be tantalizing signature of a 
discovery, the lack of it can not be easily used as an 
argument against the hypothesis. The cross-correlation 
signal is indeed much more sensitive than the autocorre- 
lation one to magnetic fields deflection and, importantly, 



to unknown experimental systematic effects. For the case 
at hand, the main responsables of the displacement up 
to 50 degrees between clusters in the data and overdensi- 
ties in the LSS are the hot-spots from the SUGAR data, 
which is the experiment among the ones considered which 
mostly suffers for a poor angular resolution, beside not 
well-understood systematics in the energy scale deter- 
mination. Indeed, limiting the analysis to the northern 
hemisphere experiments only, the data show quite a good 
cross-correlation with the local over-density of matter, 
especially within z ~ 0.02. More noticeably, they fall 
relatively close to the so-called Super-Galactic plane, as 
can be appreciated also from a comparison between Fig. [I] 
and Fig. This correspondence was already notice d by 
the authors of ref. [2l| and further assessed in ref. [34| . 
So, while the features in the cross-correlation function 
vary quite a bit excluding e.g. one dataset, the auto- 
correlation ones do not, as discussed more extensively 
in 1. 

Aware of this caveat, it is still worth exploring the con- 
sequences of assuming that displacements up to 50° with 
respect to the true sources are effective. Under this hy- 
pothesis, and, if the signal corresponds to extragalactic 
structures, we would be brought to conclude that: (i) if 
UHECRs arc dominated by protons, then there are sig- 
nificant deflections by extragalactic magnetic fields. In- 
deed, although the GMF may be not well reproduced by 
current models, even changing within reasonable ranges 
the GMF geometry and intensity no appreciable cross- 
correlation at small angles appears, (ii) If there is a 
significant fraction of heavy nuclei in the UHECR flux, 
results may be also explained with a negligible role of 
extragalactic magnetic fields, attributing to GMF deflec- 
tions the significant (~ 30°-50°) displacement between 
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FIG. 6: Chance probability of autocorrelation of the data for different GMF models and for the PH model maps (left) and SH 
model maps (right) taken as reference sample. The cases Z = 1 and Z = 2 are shown in each panel. 



the observed clusters in the data and the real galaxy 
clusters. Note that, if in any case overall deflections 
as large as ~ 50° are effective, peculiar manifestations 
of regular deflections in the magnetic field — like elonga- 
tions directed towards structures with a proper ordering 
of energies — may not be observable due to non-negligible 
non-linearities, especially in the case of a chemically in- 
homogeneous sample of UHECRs. Yet, if this interpre- 
tation would turn out to be the correct one, we expect 
that once higher statistics will be available at higher en- 
ergies, these features will eventually show up in the data, 
a prediction that can eventually be confirmed by Auger. 



IV. DISCUSSION AND CONCLUSION 



The anisotropy pattern of the combined UHECRs data 
(re-scaled energy E = 4 x 10 19 eV in HiRes scale), al- 
though compatible with isotropy at very large angular 
scales, shows a peculiar medium scale clustering corre- 
sponding to 6-7 spots of roughly 20°-30° degrees of exten- 
sion distributed uniformly in the sky. If confirmed, this 
would have a wealth of consequences for the long-awaited 
astronomy of UHECRs. At a general level, the absence of 
a correlation with the Galactic Plane or the absence of an 
excess toward the Galactic Center disfavor respectively 
Galactic astrophysical sources and heavy relic decays in 
our Galactic Halo as origin of these events. Extragalactic 
sources of astrophysical nature appear instead the most 
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FIG. 7: Cross-correlation chance probability between the dif- 
ferent LSS models and the data, with the solid lines represent- 
ing the result uncorrected for the GMF. Upper set of curves 
in each panel refer to the PH model while lower ones to the 
SH model. 



likely accelerators consistent with these features. Should 
forthcoming data show similar features, the first realis- 
tic quantities one could extract from the clustering are 
constrains the number density of the sources as well as 
their type, from their bias with respect to LSS (see [H[ 
and references therein). As next goal, one should be able 
to establish explicit cross-correlation with extragalactic 
structures and/or hints in that direction, as elongations 
of events directed toward the sources, whose distance 
should scale inversely to the UHECR energy. More "con- 
ventional" astronomy, determining the locations of single 
UHECR sources and perhaps of their spectrum is likely 
demanded to a subsequent phase when sufficient statis- 



tics will be available. 

Although the clustering properties in the data are in- 
triguing, the interpretation of the signature is puzzling, 
especially in absence of a significant statistics at higher 
energy and of chemical composition constraints. The 
comparison between significant sets of data at different 
energy cuts may reveal the importance of magnetic de- 
flection effects. At the moment, we can only speculate on 
the possible implications of the signal — assuming it is not 
a statistical fluke — under some simplifying hypotheses. 

One possibility is that these excesses may trace LSS 
overdensities in the near universe (within the GZK- 
sphere). The autocorrelation analyses reported in this 
paper show that this interpretation is indeed favored in 
particular if the effective horizon is smaller than the GZK 
one for protons of the assumed energy. Both a significant 
fraction of heavier nuclei and a significant role of extra- 
galatic magnetic fields may cause this effect (the former 
might be favored by recent Auger data [36|].) Although 
not statistically significant, this interpretation may be 
supported by a weak hint of a broad minimum in the 
cross correlation function (at the level of nominal chance 
probability of 10%-15%) around 50° if a small horizon 
(z < 0.02) is assumed. Both signatures are relatively ro- 
bust with respect to deflections in typical GMF models, 
although some marginal improvement or worsening may 
arise for some choices of the GMF model and effective 
rigidities. In this case, the size of the hot-spots would 
be due partly to the one of the largest overdensities in 
the local LSS and partly to magnetic smearing needed 
to explain the overall deflection with respect to the LSS. 
The latter effect would be in general subleading but for 
the SUGAR hotspots in the Southern Emisphere, which 
are the most distant ones from overdensities. This may 
be physically associated to the more intense magnetic 
fields towards the central regions of our Galaxy to which 
SUGAR is pointing. An alternative interpretation of the 
data is that they are due to very few (0(5-6)) powerful 
sources. Yet, the smearing of a point-like emission to 
the level of the observed spots of O(20°) would require a 
quite extreme magnetized environment [2^ . |30| . 

In any case, the hints for some structures in the data 
are very exciting, and we urge an independent cross-check 
with the nowadays large statistics collected by Auger. If 
confirmed, together with the indication for the presence 
of a GZK-like feature in the energy spectrum of HiRcs 
data this likely implies that UHECR are dominated 
by astrophysical sources (as opposed to exotic scenarios) . 
However, far from being the end of the UHECR saga, the 
combined use of spectral information, chemical composi- 
tion constraints, and anisotropy maps at different ener- 
gies would offer the tools for the long-awaited hunt for 
the UHECR accelerators. 
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